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Abstract 

The CUORE {Cryogenic Underground Observatory for Rare Events) experiment 
projects to construct and operate an array of 1000 cryogenic thermal detectors of 
Te02, of a mass of 760 g each, to investigate rare events physics, in particular, double 
beta decay and non baryonic particle dark matter. A first step towards CUORE 
is CUORICINO, an array of 56 of such bolometers, currently being installed in 
the Gran Sasso Laboratory. In this paper we report the physics potential of both 
stages of the experiment regarding neutrinoless double beta decay of ^^'^Te, WIMP 
searches and solar axions. 

PACS: 23.40.-s; 95.35.+d; 14.80.Mz 
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1 Introduction: The CUORE project 

Rare event Physics at the low energy frontier is playing a significant role in Particle 
Physics and Cosmology. Examples of such rare phenomena could be the detection of non- 
baryonic particle dark matter (axions or WIMPs), supposedly filling a substantial part of 
the galactic haloes, or the neutrinoless double beta decay. These rare signals, if detected, 
would be important evidences of a new physics beyond the Standard Model of Parti- 
cle Physics, and would have far-reaching consequences in Cosmology. The experimental 
achievements accomplished during the last decade in the field of ultra-low background 
detectors have lead to sensitivities capable to search for such rare events. Dark matter 
detection experiments have largely benefited from the techniques developed for double 
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beta decay searches. Both types of investigation, which have very close relation from an 
experimental point of view, are the two main scientific objectives of CUORE. 

Due to the low probability of both types of events, the essential requirement of these 
experiments is to achieve an extremely low radioactive background. For that purpose the 
use of radiopure detector components and shieldings, the instrumentation of mechanisms 
of background identification, the operation in an ultra- low background environment, in 
summary, the use of the state-of-the-art of low background techniques is mandatory. In 
some of these phenomena, like the case of the interaction of a particle dark matter with 
ordinary matter, a very small amount of energy is deposited, and the sensitivity needed 
to detect events within such range of energies, relies in how low the energy threshold of 
detection is. In addition, to increase the chances of observing such rare events a large 
amount of detector mass is in general advisable, and so most of the experiments devoted to 
this type of searches are planning to enlarge their detector mass, while keeping the other 
experimental parameters (background, energy thresholds and resolutions) optimized. 

On the other hand, the recent development of cryogenic particle detection |I| has lead 
to the extended use of thermal detectors to take advantage of the low energy threshold 
and good energy resolution that are theoretically expected for the thermal signals. This 
detection technology has also the bonus of enlarging the choice of materials which can be 
used, either as DM targets or 2j3 decay emitters. After a long period of R&D to master 
the techniques used in cryogenic particle detectors, low temperature devices of various 
types are now applied to the detection of double beta decay or particle dark matter 0. 
A major exponent of this development is the MIBETA {Milano Double Beta) experiment 
, which has successfully operated a 20 Te02 crystal array of thermal detectors of a total 
mass of 6.8 kg, the largest cryogenic mass operated up to now. With the objective of going 
to larger detector masses and of improving the sensitivity achieved in the smaller arrays, 
the CUORE {Cryogenic Underground Observatory for Rare Events) project was born 
some years ago as a substantial extension of MIBETA. The objective of CUORE is to 
construct an array of 1000 bolo metric detectors with cubic crystal absorbers of tellurite 
of 5 cm side and of about 760 g of mass each. The crystals will be arranged in a cubic 
compact structure and the experiment will be installed in the Gran Sasso Underground 
Laboratory. The material to be employed first is Te02, because one of the main goals 
of CUORE is to investigate the double beta decay of ^^°Te, although other absorbers 
could also be used to selectively study several types of rare event phenomenology. Apart 
from the MIBETA experiment, already completed, a wide R&D program is under way 
in the framework of CUORICINO, a smaller and intermediate stage of CUORE, which 
consists of 56 of the above crystals (total mass of 42.5 kg). CUORICINO is currently 
being mounted in the Gran Sasso Laboratory (LNGS) and is, by far, the largest cryogenic 
detector on the stage. The preliminary tests are encouraging 0. 

In the present work the prospects of CUORE and CUORICINO experiments with 
respect to their double beta decay discovery potential and their detection capability of 
WIMP and axions are presented. 
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2 Background considerations 



From the radioactive background point of view, thermal detectors are very different than, 
say, conventional ionization detectors, the former being sensitive over the whole volume, 
implying that surface impurities may play an important role. In principle, bolometric de- 
tectors like those of the MIBETA, CUORICINO and CUORE experiments are expected 
to have radioactive backgrounds larger than that of the conventional germanium ioniza- 
tion detectors, also because the more complex technology of cryodetectors had not been 
yet fully optimized from the point of view of the radiopurity of the components near the 
detector. On the other hand, the production of radiopure absorbers for thermal detectors 
is less mastered than the production of radiopure Ge crystal. However, after a consid- 
erable R&D effort a significant improvement in the radiopurity of bolometers has been 
accomplished during the last years, and now these devices have achieved very competitive 
backgrounds [0, |, g |ig]. 

As one of the main ingredients in the prospective physics potential of the CUORI- 
CINO and CUORE experiments is the level of background achievable, an evaluation of 
the expected background must be done, including the contribution of the possible back- 
ground sources. A preliminary estimation was made in the Status Report of CUORICINO 
presented to the LNGS Scientific Committee [|lT|. A more complete estimate is underway. 
In the following we will give a briefing of that (unpublished) reference and will add further 
considerations. 

There are several background sources to be considered: the environmental backgrounds 
of the underground site (neutrons from the rocks, ambient gamma flux...), the intrinsic 
radioactivities of the detectors, components and shielding (bulk and surface), the cos- 
mic muon induced backgrounds (neutrons, muon direct interactions . . . ), as well as the 
possible cosmogenic induced activities produced when the detector and components were 
outside the underground laboratory. Even small, the leakage of double beta counts emit- 
ted by the absorber's nuclei in the relevant region of analysis might be a potential source 
of background. The depth at which the experiment will be performed, plus the addition 
of an efficient cosmic veto and a suitable passive shielding (lead, polyethylene, . . . ) will 
effectively reduce the external background. So we will refer mainly to the intrinsic back- 
ground, although we will make some remarks concerning the other sources. On the other 
hand, the experience gained in the MIBETA experiment, where the measured background 
is used as a test-bench for checking the MC estimates, has been essential to know which 
could be the expectations for CUORICINO. Using the MIBETA results and the CUORI- 
CINO tests a Monte Carlo study of the expected intrinsic background has been carried 
out for the CUORICINO (and CUORE) set-ups. Supposed radioactivity impurities in the 
bulk of the tellurite crystals, as well as in the dilution refrigerator and surrounding shield- 
ing were taken into account. Background suppression due to the anticoincidence between 
crystals (which will be significant in the CUORE set up) was worked out. The effect of the 
a surface contaminations observed in the MIBETA crystals was also evaluated and scaled 
up to CUORICINO. The simulation was performed with the GEANT code. Alpha, beta 
and gamma emissions (supposed to be uniformly distributed and isotropically emitted) 
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from the natural chains {^^^U, '^^'^Th) and from other isotopes {"^^K, ^^Co, "^^^Pb) were 
included. 

In the case of MIBETA, the measurements based on the a peaks of the background 
have shown that the bulk intrinsic contamination of U/Th inside the crystals are very low 
(less than a few 10~^° g/g) and that the main sources of background are located on the 
surfaces. According to the method used to clean the holder of the crystals and to lap the 
surfaces (and depending on the radiopurity of the powders used in the lapping), the back- 
ground values, in the best of the cases, stand about roughly half-a-count per keV, kg and 
year in the 2500-2556 keV (2/3 decay) region, more precisely, b=0.68±0.06 c/keV/kg/y in 
the "old" and b=0.35±0.14 c/keV/kg/y in the "new" configuration (labelled respectively 
after the best lapping operation and a new mounting system) §]. A background spectrum 



of MIBETA corresponding to the 2/3 decay region is shown in Fig. 2 of Ref. The 
addition of a neutron shielding of 10 cm of borated polyethylene does not reduce the 
background, as expected if the dominant contribution comes from the surface. With the 
above background value the lower limit of the neutrinoless half-life of ^^°Te obtained is 
— 1-4 X 10^^ y, or equivalently, an upper bound for the Majorana neutrino mass of 
{rriu) <2 eV (when using the Heidelberg nuclear matrix element calculation of Ref. p9[]), 
which is the second-best published result |jlO|]. Similarly, as far as the low energy region 
is concerned, the most recent current background of MIBETA, obtained after the new 
lapping of the crystal detector surfaces, stands around 1 c/keV/kg/day between 10-50 
keV and 0.2 c/keV/kg/day between 50-80 keV, region where the dark matter signal is 
expected 0. Ref. shows the MIBETA background spectrum corresponding to the low 
energy, dark matter region. With these background values, the exclusion plot of WIMPs 
interacting coherently with Te and O is depicted as the (dashed) contour of figure |l]. Both 
results of MIBETA, for the low (<100 keV) and the high (~2500 keV) energy regions will 
serve as references for the CUORICINO/CUORE projections. 

To reproduce and understand the contribution to the MIBETA background, a MC 
calculation has been made for MIBETA (and then extended to the CUORICINO and 
CUORE geometries) including the main components and geometry of the experimental 
set-up (array of crystals, holder, cryostat, shielding, . . . ). The radioactive impurities used 
as inputs were all the contaminants identified in the MIBETA background spectrum, 
through the alpha or gamma lines. They are ^^^U in secular equilibrium, ^^^Bi, ^^°Pb, 
^^^Th in secular equilibrium, "^^K and ^°Co. The amount in which each impurity is present 
is obtained as the upper limit for which the MC simulation matches best the measured 
background. A previous hint about the amount of such radioimpurities was obtained by 
measuring in low background HPGe detectors (at Gran Sasso and Canfranc underground 
Laboratories) the materials and components typically used in the experimental set-up. 

The main results -which should be taken as the expectations for the CUORICINO 
background- are the following: in the very low energy region (from 10 to 100 keV), the 
background is dominated by ^^''Pb surface contamination located either on the crystal 
surface or on the surface of the copper mounting box. This contamination has been 
identified through the 46 keV peak and the 5.3 and 5.4 MeV alpha lines. The continuum 
between 100 and 1500 keV seems to be determined by a surface contamination of the 
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crystals/copper mounting box. The continuum counting rate in the energy region 1500- 
4000 keV could be explained by a crystal surface contamination of ^^^U, ^^^Ra, ^^^Th 
and ^^°Pb. This surface contamination seems to be the dominant contribution to the 
background in the 2(3 decay region, and only a minor contribution should be attributed 
to the ^'^^Th bulk contamination. 

The U/Th and ^^°Pb surface contaminations of the old MIBETA crystals (which have 
been of big concern in the R&D of CUORICINO) have been significantly reduced (but 
not eliminated) by means of a new lapping procedure of the crystals. These surface 
contaminations were estimated by the MC simulation of the MIBETA set-up and properly 
checked with the experimental data. They produced an overall rate of about ~0. 1-0.2 
c/keV/kg/y in the 2.5 MeV region and 0.2-0.3 c/keV/kg/d in the 10-100 keV energy 
region, dominating in both cases over the bulk intrinsic background. The new cleaning 
procedure carried out in the 3x3x6 cm^ crystals of MIBETA has reduced these values 
by more than one order of magnitude. 

In fact, in the CUORICINO MC simulation the following inputs have been assumed: 
the bulk contamination is the upper limit deduced for such radioimpurities in the cor- 
responding component of MIBETA. As far as the surface radioimpurities, they are as 
follows: the is assumed to be the same as in MIBETA. The U/Th surface contamina- 
tion has been scaled from that of MIBETA, but taking into account that they are factors 
1.5 and 2 smaller per unit surface than in MIBETA (as deduced from CUORICINO crys- 
tals tests). As far as the ^^°Pb surface contamination in the crystals and/or the copper, 
similar values that those derived from MIBETA have been assumed. Notice however that 
a reduction of one order of magnitude (or more) in this background is expected to be 
obtained and that very likely possibility will be taken into account to derive the physics 
prospects of the experiment. 

With the above proviso, the MC simulated contribution to the CUORICINO intrinsic 
background (bulk and surface) are b2/3(bulk)~2x 10~^ c/keV/kg/y and b2/3(surface)~10~^ 
c/keV/kg/y (from U/Th in the crystals and from ^^°Pb in the crystals and copper) 
in the 2.5 MeV region and b£)jv/(bulk)~3xl0^^ c/keV/kg/d and b^A/ (surface) ~5x 10~^ 
c/keV/kg/d (from U/Th in crystals) and 2xl0~^ c/keV/kg/d (from ^^°Pb in crystals and 
copper) in the 10-100 keV region. 

Taking into account that the surface contaminations could be reduced a factor 10, one 
can expect that the global intrinsic background of CUORICINO could be in the most 
simple and conservative estimate of ~a few 10~^ c/keV/kg/y around 2.5 MeV and a few 
10~^ c/keV/kg/d in the DM low energy region. 

Concerning the cosmogenic activation of the crystals produced by cosmic rays when 
they were above ground (during fabrication and transportation of the crystals from the 
factory to the underground laboratory) we have carried out a Monte Carlo estimation 
by using the code COSMOPACK Q. According with the history of the MIBETA (and 
CUORICINO) detectors, we have assumed 2 months of exposure to cosmic rays and 
then 1 year of cooling down of the induced activities. The radionuclei produced by the 

^We thank Y. Ramachers for providing us the code used by the CRESST CoUaboration 



5 



activation of tellurium are mostly tellurium isotopes (A=121,123,125,127) as well as ^^^Sb, 
^^^Sb and tritium, these last two being of more concern because of their half-life (beta 
decays of 2.7 years, end-point energy of 767 keV and 12.3 year, end-point energy of 18 
keV, respectively). The induced ^^"^Sb has also a beta decay (half-life of 60 days, end-point 
energy of 2905 keV) and a few gamma lines. The total induced activities remaining after 
one year underground have been estimated to be 6.5 c/kg/d (^^^Sb), 1.8 c/kg/d (^^^Sb) 
and 0.1 c/kg/d ('^H). In the case of tritium, for instance, the contribution to the event 
rate (from threshold (5 keV) onwards) is 5 x 10"'^ c/keV/kg/d. Their incidence in the 
background count rate are still below the intrinsic radioactivities (bulk and surface) of 
the crystals but they might set a bound for the achievable background of the experiment, 
unless the cooling down period be extended. No special incidence of the cosmogenic 
activation has been seen in the recent MIBETA spectra. However, some lines of the 
above mentioned isotopes have been observed when the crystals were temporarily taken 
out of LNGS for lapping the surfaces and operated again. 

Notice that, in any case, the cosmogenic activation is already included in the ex- 
perimental data of the MIBETA background and so, it has been taken into account in 
the extrapolation made in going from MIBETA to CUORICINO (both types of crys- 
tals follow, roughly, the same history of fabrication in China, transportation and storage 
underground) . 

On the other hand, the leakage of 2f3 counts in the low energy region (< 100 keV), 
due to the double beta decay of ^^°Te is totally negligible (< 10~^ c/kev/kg/d, assuming 
for the ^^°Te a 2f]2u half-life of ~ 10^^ y, an average of the two current measured results 
of 7.3x1020 y and 1.2x10^1 y 

As noticed before, neither cosmic muons nor neutrons have been taken into account, in 
detail, in the estimation of the background. However, the following simplified arguments 
will serve to have an approximate idea of which could be their contribution. 

The depth of the LNGS (3500 m.w.e) reduces the muon flux down to ~ 2x10^^ 
cm~2s~^, but a further effective reduction can be obtained with the use of an efficient 
(99.9%) active veto for not to miss muons traversing the detectors and to tag possible 
events associated with them. Consequently, the muon-induced events contributing to the 
background are expected to be only a small component of it. The main contribution to the 
expected background, as seen in MIBETA, is coming from the intrinsic radioimpurities 
(bulk or surface). 

On the other hand, the shielding will substantially reduce the event rate due to par- 
ticles external to the detector from various sources (neutrons and photons), from ra- 
dioactivity in the environment (natural decay chains U/Th, ^^^Pb, ^°K, . . . ), as well as 
muon-induced in the surroundings or in the shielding itself. The passive shielding typi- 
cally consists of a neutron screen (blocks of -borated- polyethylene of 10-20 cm thickness) 
to attenuate and moderate neutrons and a shell in lead (of about 10-15 cm) to attenuate 
the incoming external photons; the innermost part is made of archaeological lead (2000 
years old) with a very small content (<4 mBq/kg) of ^^^Pb (half-life 22 years) and of 
high radiopurity. The successive shielding barriers, active and passive, guarantee a very 
substantial reduction of the external background sources (we will show how this reduction 
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applies to the CUORICINO case). 

Neutrons may constitute a worrisome background in dark matter experiments because 
for appropriate neutron energies (few MeV) they can produce nuclear recoils (^100 keV) in 
the detector target nuclei which would mimic WIMP interactions. Simple kinematics tells 
that in the case of tellurium, neutrons of 1(5) MeV could elastically scatter off tellurium 
nuclei producing recoils of energies up to 31 (154) keV. In general, one considers neutrons 
of two origins: from radioactivity in the surroundings or muon-induced. Depending on 
the overburden of the underground site (i.e., depending on the muon flux), muon-induced 
neutrons are produced, at lesser or greater rate, both inside and outside the shielding. 
They are moderated (according to their energies) by the polyethylene /lead shield (when 
produced outside) or tagged by the muon veto coincidence (when produced within the 
passive shielding). A certain fraction of the neutrons of, say, a few MeVs produced outside 
can pass the veto reaching the detector and producing "dangerous" nuclear recoils and 7 
background. However, a significant fraction of the associated events can be vetoed because 
of their interaction with the veto and also because of the hadronic showers trained by the 
muons (see Ref. 0). We have MC-estimated -see later- how many of these neutrons can 
punch through the shielding. 

In the case of external neutrons (from the rocks, from fission processes or from (n,a) 
reactions, as well as neutrons originated by muons in the walls of the underground site), 
the ambient neutron flux has been measured in LNGS. The result is of ~ 1 x 10~^ cm~^s~^ 
for the thermal component, ~ 2 x 10~® cm~^s~^ for the epithermal and ~ 2 x 10"'' cm~^s~^ 
for energies over 2.5 MeV [^]. They are fairly well moderated by the polyethylene and 
eventually absorbed or captured. We have carried out a Monte Carlo simulation of the 
propagation of neutrons through a typical (but simplified) shielding of CUORICINO. In 
particular, we have estimated the distance distributions for neutrons (of energies 1 and 
5 MeV) to thermalize (E<0.1 eV) in polyethylene. The maxima of the distributions are 
at 5 and 7 cm respectively. For instance, 90% of neutrons of 1 MeV (5 MeV) thermalize 
after 12 cm (22 cm) of polyethylene. In a shielding of 40 cm of polyethylene, 99.7 % 
of neutrons of 5 MeV are moderated down to 0.1 eV (as well as practically all neutrons 
of 1 MeV). For neutrons of higher energies the fraction of them which are thermalized 
after, say, 40 cm of water is ~92% (for neutrons of 10 MeV), ~83% (for neutrons of 25 
MeV) and ~50% (for neutrons of 50 MeV). Then, taking into account the energies of 
fission neutrons (rocks, . . . ) and the small number of cosmic muon-induced neutrons at 
3500 m.w.e., a neutron shielding of ~ 40 cm of polyethylene (or water) would practically 
reduce the external neutron flux down to energies not dangerous from the point of view 
of the nuclear recoils they could produce. 

On the other hand, the incoming flux of external neutrons becomes attenuated in 
the shielding material (in the polyethylene neutron moderator or in the gamma lead 
shielding). As said before, we have developed a MC simulation of such attenuation as well 
as the associated photon production originated by these neutrons (via inelastic scattering, 
radiative capture) in the shielding material. The MC estimate assumes a shielding box of 
external dimensions typical of those of CUORE (lmxlmx2m), with walls of lead 
and polyethylene of various thickness, to determine the fraction of incident neutrons (of 
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initial energies of 1 and 5 MeV) which are able to punch through the thickness of the 
shielding box. Taking, for instance, a box consisting of 10 cm of polyethylene (external) 
and 10 cm of lead (internal), it turns out that about 2%, 4% and 10 % of external 
neutrons having respectively energies of 0.1, 1 and 5 MeV, are able to punch through that 
shielding. The number of neutron-induced photons, which arrive inside the box (in the 
detector volume) has been also evaluated. Starting, for instance, with neutrons of 5 MeV 
entering the shielding and reaching the region of the crystals with energies of 2.6 MeV 
(generated by inelastic scattering with ^°^Pb) or with 2.2 MeV (due to thermal neutron 
capture in hydrogen), they produce 10~^ photons per neutron (or 8 x 10^^ photons per 
neutron) respectively. In the case of neutrons of 1 MeV, these figures are, respectively, 
10""^ photons per neutron and practically zero photons. These considerations referred to 
the typical shielding of the CUORICINO/CUORE experiment give an idea of the limited 
incidence of the external neutrons, as previously stated. 

Other neutrons, produced by muon interaction inside the shielding materials, are very 
scarce and tagged as events coincident with the muon veto. As it is well-known, muon- 
induced neutrons are originated in a variety of processes. The reduction of muon flux in 
underground sites results in a substantial depletion of the associated neutrons and below 
~100 m.w.e. the dominant sources of neutrons are nuclear fission processes and (n,a) 
reactions in rocks and other environmental material with sizeable content of U/Th. The 
energy spectrum of muon-induced neutrons is approximated by an inverse energy power 
law (^E~^'^^ for 1-50 keV and above 50 keV) but other neutron energy spectra have 
been proposed. The neutron yield per muon can be approximately evaluated through the 
simple expression Nn = 4.14 X X 10"^ neutrons/ (gcm-2) per muon which fits 

the value of 1.5x10"^ neutrons/ (gcm^^) of the muon-induced neutron flux measured by 



the LVD experiment at Gran Sasso ||T8|. Other authors have used similar expressions. 

Using an average neutron yield per muon of 2.2x 10"'^ neutrons/ (gcm~^), together with 
the LNGS muon flux (2.5x10"*^ yu/(cm^s)), that would produce in the simplified CUORE 
shielding referred above (box of 1 x 1 x 2 m^, with walls of polyethylene (10 cm) and 
lead (10 cm)), about ~0.5 neutrons/day in the polyethylene shield and ~4 neutrons/day 
in the lead shell. So, independently of the mechanism set up to reject or tag the events 
associated to neutrons, their rather small number is expected to play a secondary role 
in the total background compared with other, intrinsic, sources of background. On the 
other hand, much in the same way, the remarks we made above about the cosmogenic 
induced activity when going from MIBETA to CUORICINO also apply to the neutron 
induced background, which is already incorporated to the MIBETA background data, and 
so taken into account in the extrapolation. A more complete and detailed study of the 
external, non intrinsic background rates is underway. 

The 7-ambient background in Gran Sasso is ~ 1 cm~^s^^ [l^ and it can be attenuated 
by a proper lead shielding. 

The background values quoted for MIBETA were obtained nevertheless with a more 
simplified shielding (no active veto, no neutron shielding except in the last running, where 
10 cm of polyethylene was added). In spite of this fact, the MIBETA background 0] at 
low energy (~1 c/keV/kg/d at threshold -10 keV- and at 3500 m.w.e.) is similar than 
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the measured event rate anticoincident with the veto in the CDMS experiment |0| (~ 2 
c/keV/kg/d at 10 keV and at 20 m.w.e.) and in the EDELWEISS experiment 1.8 
c/keV/kg/d at 30 keV and 4000 m.w.e.), without veto, (in both cases, obviously, prior to 
charge-heat discrimination) and is roughly equal to that of DAMA and ANAIS |13| 
(~ 1.5 c/keV/kg/d, at 2 keV and 3500 m.w.e. and at 4 keV and 2450 m.w.e. respectively) 
but still one order of magnitude worse than that of IGEX |]14| (~ 0.2-0.05 c/keV/kg/d 
at 4 keV and 10 keV respectively, and at 2450 m.w.e.). In the double beta decay region 
(~ 2.5 MeV) the MIBETA background values (~ 0.3-0.6 c/keV/kg/y in, respectively, the 
old and new set-up) are competitive but still higher than those of the IGEX experiment 
[p!5| (0.05 c/keV/kg/y at 2 MeV) (which uses Pulse Shape Discrimination). 

The challenge of CUORICINO (and later on, of CUORE) is to reduce significantly the 
MIBETA background values (both in the low and high energy regions) by, say, two orders 
of magnitude without using background discrimination mechanisms like the simultaneous 
measurement of charge (or light) and heat. We argue that it can be achieved in two steps. 

Starting from the current background values of MIBETA (~ 1 c/keV/kg/d from 10-50 
keV, ~ 0.2 c/keV/kg/d between 50-80 keV and ~ 0.3 c/keV/kg/y around 2.5 MeV) and 
taking into account the preceding discussion about the incidence of the various compo- 
nents of the background, we expect that these background values can be reduced by one to 
two orders of magnitude, going down to values smaller than 0.1-0.05 c/keV/kg/day in the 
low energy region and to about 0.01 c/keV/kg/y (or lower) in the double beta region, as 
suggested by the MC simulation and the nature and location of the background radioac- 
tivities. That extrapolation is very conservative if one takes into account the exhaustive 
radiopurity selection of the detector materials and components which is being done for 
CUORICINO, the new lapping procedure of the crystal surface, the closer packing of the 
crystals, the suppression of significant amount of holding material, the powerful anticoin- 
cidence rejection provided by the segmented geometry of the crystal array and the use of 
a much better shielding than that employed in MIBETA. Nevertheless, as a first step in 
the background achievements, we will assume in the following prospective analysis of the 
CUORICINO capabilities, even more conservative background values (including intrinsic 
and external), more precisely of 1 c/keV/kg/day and of 0.1 c/keV/kg/day in the low en- 
ergy region (at threshold, ~5 keV) and values of 0.1 c/keV/kg/y and of 0.01 c/keV/kg/y 
in the 2-2.5 MeV region. The forthcoming results of the CUORICINO experiment will 
tell us how much this reduction can be further pursued. Even at these conservative values, 
CUORICINO is, as we will proof in the following, a powerful instrument to look for 2(3 
decays and WIMPs. In the case of CUORE the background values that will be assumed 
(also approximatively) are 0.1 and 0.01 c/keV/kg/day in the low energy region and 0.1 
and 0.01 c/keV/kg/y around 2.5 MeV, which should be considered, also, two successive 
steps of the experiment. 

Regarding the expected threshold and resolution, in the CUORICINO tests |^ energy 
thresholds of ~ 5 keV have been obtained and a resolution of 1 keV at the 46 keV line 
of ^^°Pb achieved, in some of the detectors. Ref. |0 illustrates this point. In the case 
of CUORE thresholds of 5 keV and energy resolutions of 1 keV at low energies will be 
assumed. As far as the energy resolutions obtained in the double beta decay region, values 
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of 3 keV at 2615 keV were achieved in some crystals but they are worse (by a factor two) 
in others. See Fig. 1 of Ref. |ll]]. Nevertheless, they are clearly better than that obtained 
in the previous 6x3x3 cm^ (MIBETA) crystals (8~10 keV). Taking into account these 
expectations, we discuss in the following the prospects of CUORICINO and CUORE for 
double beta decay searches (section 2), for WIMP detection (section 3) and for solar axion 
exploration (section 4). 



3 Double beta decay 

One of the main scientific objectives of the CUORE detectors is to search for the double 
beta decay of the ^^°Te isotope contained in the (natural) Te02 set of crystals. 

The importance of the nuclear double beta decay as an invaluable tool to explore 
particle physics beyond the Standard Model has been repeatedly emphasized and widely 
reported In the Standard Model of Particle Physics neutrinos are strictly mass- 

less, although there is no theoretical reason for such a prejudice. On the experimental 
side, moreover, there exist strong evidences from atmospheric neutrino data (from Su- 
per Kamiokande) and from experiments with solar neutrinos from Homestake, Gran Sasso 
and Kamioka, since long time ago, which suggest that neutrinos have indeed masses and 
oscillate among the three species. The recent results of the solar u experiment from 



SNO with both CC (Charged Current) and NC (Neutral Current) interactions also 
combined with SuperK have provided, definitively, a strong evidence that neutrinos do 
oscillate and, consequently, the existence of non-zero mass neutrinos. However, neutrino 
oscillation experiments provide the squared mass difference between the neutrino species 
but not their absolute value and scale. The neutrinoless double beta decay would help to 
solve this question and to disentangle the hierarchy scheme of the neutrino flavours. Most 
of the models (see ||22|) indicate that the Majorana neutrino mass parameter could be 
around (or slightly below) (m^) ~ 0.05 eV, value within reach of the future double beta 
decay experiments, like CUORE. On the other hand, galaxy formation requires a small 
amount of hot non-baryonic dark matter likely in form of neutrinos to match properly the 
observed spectral power at all scales of the universe. The question of the neutrino mass 
is one of the main issues in Particle Physics. 

In the Standard Model, neutrinos and antineutrinos are supposed to be different parti- 
cles, but no experimental proof has been provided so far. The nuclear double beta decay 
addresses both questions: whether the neutrinos are self-conjugated and whether they 
have non-zero Majorana masses. In fact, the lepton number violating neutrinoless double 
beta decay {A, Z) {A, Z + 2) + 2e~ (2/90z/) is the most direct way to determine if 
neutrinos are Majorana particles. Moreover, the observation of a 2(30i> decay would imply 
a lower bound for the neutrino mass, i. e. at least one neutrino eigenstate has a non-zero 
mass. 

Another form of neutrinoless decay, {A, Z) {A, Z + 2) + 2e^ + x niay reveal also 
the existence of the Majoron (x), the Goldstone boson emerging from the spontaneous 
symmetry breaking of B — L, of most relevance in the generation of Majorana neutrino 
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masses and of far-reaching implications in Astrophysics and Cosmology. These and other 
issues, make the search for the neutrinoless double beta decay an invaluable tool of ex- 
ploration of non-standard model physics, probing mass scales well above those reached 
with accelerators. That is the motivation why there are underway dozens of experiments 
looking for the double beta decay of various nuclei [^] like ^^Ge (IGEX [^], Heidelberg- 



Moscow |2|), Mo and others (NEMO H, ELEGANTS H) and ^'^°Te (MIBETA, 
CUORICINO) and a few big experimental projects, like CUORE, Majorana |1 C^Ge), 



MOON 13 (i°°Mo) and EXO |28| (136Xe). 

The cryogenic thermal detectors provide new double beta emitter nuclei to be explored 
in "active" source=detector calorimeters. Some of them have been tested and others are 
already in running detectors, like '^^Ca in CaF2, ^^^Te in Te02, and ^^^Cd in CdW04. As 
far as the Tellurium Oxide is concerned, the 130-Tellurium isotope is a good candidate for 
double beta decay searches: its isotopic content in natural Tellurium is 33.87%, and its 2/? 
Q-value (Q2/3 = 2528 ±1.3 keV) is reasonably high to escape from the main radioimpurity 
lines when looking for a neutrinoless signal. Moreover, this Q-value happens to be between 
the peak and the Compton edge of the 2615 keV line of ^''^Tl, which leaves a clean window 
to look for the signal. Finally, it has a fairly good neutrinoless nuclear factor-of-merit 
= Gou \M^'^\^, where Gq^ is an integrated kinematical factor qualifying the goodness 
of the Q2/3 value (large phase space) and M^'^ the neutrinoless nuclear matrix elements 
characterizing the likeliness of the transition. 

In table |1| we quote the values of in the ^^^Te case calculated in various nuclear 
models [Q, together with those of other emitters used in source=detector calorimeters. 
It can be seen that no matter the nuclear model used to compute the neutrinoless decay 
matrix elements, the merits of ^^^Te are a factor 5-10 more favorable than those of ^^Ge 
(the emitter where the best neutrinoless double beta decay half-life limits have been 
achieved so far), which translates into a factor 2 to 3 better as far as the (m^) (Majorana 
neutrino mass parameter) bounds are concerned. 

The detector factor-of-merit F^, or detection sensitivity, introduced long time ago 
by Fiorini, provides an approximate estimate of the neutrinoless half-life limit achievable 
with a given detector. For source=detector devices, it reads: 



F°^ = 4.17x10^^x^^^X6 years (1) 

where A is the atomic mass, a is the isotopic abundance, M the detector mass in kg, b 
the background in c/keV/kg/y in the 2/3 neutrinoless decay region, t the running time 
in years, F the FWHM energy resolution in keV and e the detector efficiency (which is 
practically one within the fiducial volume of the detector). In the case of a Te02 crystal 

detector, F^^ ~ 8.86 x 10^^ -y/^^, with M the crystal mass in kg and b the background in 
counts per keV and year per kg of detector mass. 

The simplest projection for CUORICINO uses a background of b=0.3 c/keV/kg/y (at 
2.5 MeV) and an energy resolution FWHM(2.5 MeV)~ 8 keV, which are the performances 
obtained currently in MIBETA (notice that the background in the "old" MIBETA ver- 
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sion was of ~0.6 c/keV/kg/y). With that proviso, the sensitivity of 1 year exposure of 
CUORICINO (42.5 kg mass of crystals) will be T^J^ > 3.7 x lO^^ y. As far as the (m^) 

limit is concerned (given as usually by (m^) < mj ^/F^Ff ), since the nuclear factor-of- 
merit of ^^°Te is larger than that of ''^Ge (see Table |l|) the ^^°Te emitter provides a factor 
2-3 better {mj) bound than that obtained from ^^Ge for the same half-life limit. In other 
words, one needs to reach a half-life limit of only T^j2 > 5 x 10^^ years to get the same 
{rrii,) upper bound currently achieved in the best Germanium experiments (1.6 — 1.9 x 10^^ 
years respectively and [^). In terms of {rrty) bound, the above simplest projection of 
T^j2 — 3-7 X 10^"^ y means (m^) < 0.36 — 0.21 eV using |2^ or [^, i.e., values comparable 



to, or better than, the best currently achieved with Ge experiments. 

A second, still very conservative projection would be provided in the case of reaching 
a background of b=0.1 c/keV/kg/y (as discussed in Section 2) and a resolution of F = 
5 keV in the 2.5 MeV region. Such performances can be easily achieved in the very first 
stage of CUORICINO (notice that MIBETA has already achieved 0.3 c/keV/kg/y and 
that the energy resolution in the tests of the CUORICINO crystals ranges from F = 3 keV 
to ~ 6 keV at 2615 keV, i.e., notably better than that of smaller 3x3x6 cm^ of the 20 
crystal array of MIBETA). In that case, one would have ~ 1.25 x 10^^ a/ Mt years, with 
M (kg) the mass of the Tellurite crystal array. For the mass of CUORICINO (M=42.65 
kg), one will have a sensitivity of Fd ~ 8.15 x 10^^\/t years. Using a typical average 
value of Fiv = 4 X lO^^^^ y-\ as obtained in QRPA |g, |31|, |3|, |5[, CUORICINO 



will have a mass bound sensitivity of {rriiy) < 0.28 eV in one year, in the least favourable 
case. Using, in particular, (for comparison purposes) the value of Fn = 5.33 x 10~^^ y~^ 
of Ref. which is usually employed in the (m^) bound (~0.33 eV) derived from the 
Ge experiments, CUORICINO would provide with these assumptions (mj,) <0.24 eV. 

To go further, one needs to increase the mass of Te02 (CUORE 760 kg) and to reach 
even a lower background, which the anticoincidence capability of the CUORE array of one 
thousand crystals and the surface radiopurity could, likely, achieve going down to, say, 
b=0.05-0.01 c/keV/kg/y. In the case of the full CUORE detector, formed by one thousand 
crystals of a total mass of 760 kg, even keeping, conservatively, b and F as above (b=0.1 
c/keV/kg/y, F(2.5 MeV)=5 keV), one would get T^J^ > 3.4 x lO'^^y/i years, which in one 
year of statistics would provide {m,y) bounds ranging from 0.12 eV 0.14 eV [0, 0.26 
eV ||3^ or 0.07 eV just to mention a few nuclear matrix element estimates. However, 
the R&D to be carried out in CUORE and the operation in the antincoincidence mode is 
expected to provide better figure-of-merit than the values used in the previous predictions. 
According to such expectatives, values of 6 ~ 0.01 c/keV/kg/y and a FWHM energy 
resolution of F ~ 1 keV will be used for the last step of CUORE, i.e., a detection sensitivity 
of Ff) ~ 8.86 X 10"^^^/ Mt years, from where one can work out the best expectatives of 
the whole CUORE (M = 760 kg of tellurite), i.e., Fd ~ 2.5 x lO^^Vt years, or (m^) 
bounds ranging from ~ 0.05 t'^/^ eV (in m, mmM,m m) to ~ 0.03 t'^/^ eV (in 



32| , |33|). The ultimate sensitivity of CUORE (as far as the Majorana neutrino mass 



bound is concerned) stagnates at ~ 0.05 eV with a very softened dependence with time. 
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4 WIMP detection 



Recent cosmological observations |39] provide compelling evidence for the existence of 
an important component of non-baryonic cold dark matter in the Universe. Among the 
candidates to compose this matter, Weakly Interacting Massive Particles (WIMPs) and 
axions are the front runners. The lightest stable particles of supersymmetric theories, like 
the neutrahno describe a particular class of WIMPs. 

Under the hypothesis of WIMPS as main component of the dark matter, these parti- 
cles should fill the galactic haloes and explain the flat rotation curves which are usually 
observed in many galaxies. The detection of such particles could be attempted both by 
means of direct and indirect methods. The direct detection of WIMPs relies on the mea- 
surement of their elastic scattering off the target nuclei of a suitable detector 0. The 
non relativistic and heavy (GeV - TeV) WIMPs could hit a detector nucleus producing a 
nuclear recoil of a few keV. Because of the small WIMP-matter interaction cross sections 
the rate is extremely low. In the case of SUSY WIMPs most of the cross section pre- 
dictions [^, (derived using MSSM as the basic frame implemented with different 
unification hypothesis) encompass a range of values several orders of magnitude wide (the 
so-called scatter plots) providing rates ranging from 1 c/kg/day down to 10^^ c/kg/day 
according to the particular SUSY model. 

It is well known that the predicted signal for the WIMP elastic scattering has an 
exponentially decaying energy dependence, hardly distinguishable from the background 
recorded in the detector. The simple comparison of the theoretical WIMP spectrum with 
the one experimentally obtained, provides an exclusion (at a given confidence level), as 
dark matter component of the halo, of those WIMPs with masses (m) and cross sections 
on nucleons (a) which yield spectra above the measured experimental rate. To claim 
a positive identification of the WIMP, however, a distinctive signature is needed. The 
only identification signals of the WIMP explored up to now are provided by the features 
of the Earth's motion with respect to the dark matter halo. In particular, the annual 
modulation is originated by the combination of the motion of the solar system in 
the galactic rest frame and the rotation of the Earth around the Sun. Due to this effect, 
the incoming WIMP's velocities in the detector rest frame change continuously during 
the year, having a maximum in summer and a minimum in winter. Therefore the total 
WIMP rate changes in time with an oscillating frequency which corresponds to an annual 
period and a maximum around the beginning of June. 

The relative variation of the signal is small (a few percent) so in order to detect it 
one needs large detector masses to increase statistics and several periods of exposure to 



minimize systematics. Several experiments have already searched for this effect [^, |4g, ^ 



and since 1997 one group has reported a positive signal [O which has been appearing 
along four yearly periods. The present situation is no doubt exciting: on one hand 
that result has triggered an intense activity in the field; on the other, the experimental 
sensitivities of various types of underground detectors are entering the supersymmetric 
parameter space |^ and in particular three of them 0, H, [l^ are excluding, to a larger or 
shorter extent, the region of mass and cross-section where the reported WIMP is supposed 
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to exist. New data from one of them [^] have excluded totally the DAMA region. We will 
discuss in the following the capabilities of CUORICINO and CUORE to exclude WIMPs 
using the total time-integrated experimental rate and comparing it with the predicted 
nuclear recoil rate. To look for the annual modulation signal in CUORICINO / CUORE 
experiments, which in principle have enough masses to be sensitive to it, one needs to 
know their stability performances. The analysis of the CUORE / CUORICINO potential 
for annual modulation searches will be performed -following statistical consideration- (see 
Ref. [0), with the proviso that systematic uncertainties are under control. Data on the 



stability of MIBETA and on the first running of CUORICINO will be crucial to asses 
such hypothesis. 

To calculate the theoretical WIMP rate, standard hypothesis and astrophysical pa- 
rameters are assumed, i.e., that the WIMPs form an isotropic, isothermal, non-rotating 
halo (the isothermal sphere model) of density p = 0.3 GeV/cm^, which has a maxwellian 
velocity distribution with Vrms = 270 km/s (with an upper cut corresponding to an escape 
velocity of 650 km/s), and a relative Earth-halo velocity of Vr = 230 km/s). Other, more 
elaborated halo models, which have been considered recently would lead to different 



results. The same applies when other astrophysical parameters are employed or when 
uncertainties in the halo WIMPs velocity distribution are included The theoretical 
predicted rate is expressed in terms of the mass and cross-section of the WIMP-matter 
interaction. The cross sections are normalized per nucleon assuming a dominant scalar in- 
teraction, as is expected, for instance, for one of the most popular dark matter candidates, 
the neutralino: 

I^W,n 

where A is the target (oxigen and tellurium) mass number, iJi%/j^ is the WIMP-nucleus 
reduced mass, and the WIMP- nucleon reduced mass. The Helm parameterization 
[pO| is used for the scalar nucleon form factor. The (m, a) exclusion plot is then derived 
by requiring the theoretically predicted signal for each m and a in each energy bin to be 
less than or equal to the (90% C.L.) upper limit of the (Poisson) recorded counts. The 
bin width is assumed to be equal to the detector resolution. 

In figure |l] the exclusion plots for coherent spin-independent WIMP-matter interac- 
tion are shown for two possible values of the background of CUORICINO, 1 and 0.1 
c/keV/kg/day. The first value is of the order of the background already achieved from 
threshold onwards (10-50 keV) in the MIBETA latest results (see 0). The value 0.1 
c/keV/kg/day is a one-order-of-magnitude extrapolation from that currently achieved in 
MIBETA (see discussion on Section 1) and is close to the one obtained above 50 keV. 
The main challenge of this hypothesis is to get b=0.1 c/keV/kg/day below 10 keV. Notice 
that these values are more conservative than that derived from the MC simulation, and 
could be taken as moderate extrapolations of the background of MIBETA and of the 
CUORICINO test crystals. In the case of CUORE, background values of 0.1 and 0.01 
c/keV/kg/day will be assumed. Notice, moreover, that values of a few 0.01 have been 



obtained above 10 keV in the raw spectra of Germanium experiments (like IGEX W% ) 
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without using mechanisms of background rejection, and so it does not seem impossible to 
achieve such equivalent (< 0.04 c/keV/kg/day) small values in crystal thermal detectors 
of Tellurium (only phonons, and no discrimination mechanism). To draw the two exclu- 
sion contours of Fig. |l], a low energy resolution of 1 keV and an energy threshold of 5 keV 
have been assumed as well as an exposure of 2 years of CUORICINO (84 kg-year). The 
projected exclusion contours are compared with the one currently obtained from MIBETA 
(dashed line). In figure ||, the exclusions for the two quoted values of the background of 
CUORE, 0.1 and 0.01 c/keV/kg/day, are similarly presented for an exposure of 1 year 
(760 kg-year). In both Figures |l] and the region corresponding to the annual modulation 
positive signal reported by the DAMA collaboration is depicted as the closed "triangular" 
contour. One can see that the DAMA region could already be explored by the CUORI- 
CINO experiment with a background of 0.1 c/keV/kg/day. Of course CUORE will be 
able to explore a larger region of the (m, a) plane entering substantially into the scatter 
plots of the theoretical predictions of the various SUSY models (see and references 
therein). The reader could make his own choice according to the credibility he might 
give to the above extrapolation of the background and other parameters from the current 
obtained performances. 

As previously noted, CUORE and to some extent CUORICINO have detector masses 
large enough to search for the annual modulation signal. As it is well known, an essential 
requirement to estimate the prospects of any detector to search for annual modulation 
is to have a superbe control of systematic errors and to assure that the stability of the 
various experimental parameters, which might mimic periodic variations of the signals, 
are kept within a small percent of the (already tiny) expected signal. The various changes 
of set-up, crystals and shielding of the MIBETA experiment have not provided a definitive 
estimation of the long-term stability parameters of MIBETA. Possible instabilities are that 
of the electronic gain and the ensuing time fluctuation of the energy scale (both in energy 
thresholds and energy resolutions), the temperature variations, the possible fluctuation 
in time of the efficiency with which the triggered noise is rejected and others. They must 
be kept well below the small expected seasonal modulation of the WIMP signal. The fact 
that we are dealing with a very small signal depending on time, which typically amounts to 
a fraction between 1% and 7% of the average count rates, reinforces the need for a control 
of the stability of the experiment (well below that range) over long periods of time. In 
the case of assuming that all these fluctuations are controlled well below the levels needed 
(<1%), then one can proceed to analyze the sensitivity of CUORICINO/CUORE to the 
annual modulation signal on purely statistical grounds. This has been first attempted 
5T[| and [^, but a more extensive and rigorous approach is followed in ref. |]53 



m 



where sensitivity plots for several types of detectors (and experimental parameters) are 
presented, and in particular, for CUORE and CUORICINO. 

The sensitivity of a given experimental device to the annual modulation signal (accord- 
ing to the detector material employed and the experimental parameters of the detectors) 
has been extensively studied in Ref. on purely statistical grounds. Following the 
guidelines of that reference, it can be precisely quantified by means of the S parameter, 
defined from the likelihood function or, equivalently, from the function of the cosine 
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projections of the data (for further details see ref. p3|): 



s = = 0) - vmin ^ x {(y = 0) - xLn- (3) 

This parameter measures the statistical significance of the modulation signal detected 
in an experimental set of data. However, for a given (m, a) and a given experiment the 
expected value (5^) can be estimated using the expression derived in ref. |p3| : 



(5^) = i V ^JIlA^iII^^f^MTa + 2. (4) 
2 ^ Ok + oo,fc 

where Sm,k and 5*0,^ are the modulated and non-modulated parts of the WIMP signal in 
the fcth energy bin of AEk width, hk is the background in that energy bin and MTa the 
effective exposure, being a a coefficient accounting for the temporal distribution of the 
exposure time around modulation maxima and minima [a = 1/n ^17=1 ^'^^'^ ^i^i ~ ^o) 
n temporal bins). 

Using this equation we have estimated the region that could be within reach for 
CUORE and CUORICINO with the above mentioned assumptions on the background 
levels. We have fixed a value of 5.6 for (5^) that corresponds to 50% probability of ob- 
taining a positive result at 90% C.L.. In figure |] it is shown the curves obtained for a 
threshold of 5 keV, two years of exposure with CUORICINO (84 kg year) and two as- 
sumed fiat backgrounds of 1 and 0.1 c/keV/kg/day. One can see that CUORICINO could 
already explore pretty well the DAMA region looking for a positive annual modulation 
signal. In figure ^ similar curves are presented, assuming fiat backgrounds of 0.1 and 0.01 
c/keV/kg/day, two years of exposure of CUORE (1500 kg year) and a threshold of 5 keV 
(solid lines). The possibility of a lower thresholds of 2 keV with a background of 0.01 
c/keV/kg/day is also shown (dashed line). 

In conclusion, CUORE and CUORICINO will be able to explore and/or exclude 
WIMPs lying in large regions of their parameter space. The capability of CUORICINO 
/ CUORE to investigate the DAMA region through the exclusion plot (time integrated 
method) relies in getting a background of 0.1 c/keV/kg/day from 5 keV onwards, in- 
dependently of more elaborated time modulation methods which require an exhaustive 
control of the stability of the experiment. However, CUORICINO and CUORE could also 
attempt to look for annual modulation of WIMP signals provided that the stability of the 
experiment be guaranteed. 



5 Solar axion detection 

Axions are light pseudoscalar particles which arise in theories in which the Peccei-Quinn 



U(l) symmetry has been introduced to solve the strong CP problem They could 

have been produced in early stages of the Universe being attractive candidates for the 
cold dark matter (and in some particular scenarios for the hot dark matter) responsible 
to 1/3 of the ingredients of a fiat universe. Dark matter axions can exist in the mass 
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window 10 "^^^^ eV < < 10 ^ eV, but hadronic axions could exist with masses around 
the eV. 

Axions could also be copiously produced in the core of the stars by means of the 
Primakoff conversion of the plasma photons. In particular, a nearby and powerful source 
of stellar axions would be the Sun. The solar axion flux can be easily estimated |56| 
within the standard solar model, resulting in an axion flux of an average energy of about 
4 keV that can produce detectable X-rays when reconverted again in an electromagnetic 
field. Moreover, it has been pointed out recently that the dimming of supernovae SNIa 
might be due to the conversion of photons into axions in the extra-galactic magnetic 
field ||5^. That 7-a oscillation could make unobservable about 1/3 of the SN emitted 



light and so, they would appear fainter than implied by the luminosity-distance versus 
redshift relation, without need to invoke an accelerated expansion of the Universe. The 
SN result would be matched by axions of mass ~ 10~^^ eV and coupling to photons 
Qa-yy ~ 2.5 X 10^^^ GeV^^. So stellar axions may play an important role in Cosmology. 
We would like to stress that, although we focus on the axion because its special theoretical 
motivations, all this scenario is also valid for a generic pseudoscalar (or scalar) particle 
coupled to photons [0. Needless to say that the discovery of any type of pseudoscalar 
or scalar particle would be extremely interesting in Particle Physics. We will keep our 
discussion, however, restricted to the case of solar axions. 



Crystal detectors provide a simple mechanism for solar axion detection ||59|, ^ . Axions 



can pass in the proximity of the atomic nuclei of the crystal where the intense electric field 
can trigger their conversion into photons. The detection rate is enhanced if axions from 
the Sun coherently convert into photons when their incident angle with a given crystalline 
plane fulfills the Bragg condition. This induces a correlation of the signal with the position 
of the Sun which can be searched for in the data and allows for background subtraction. 
The potentiality of Primakoff conversion in crystals relies in the fact that it can explore 
a range of axion masses {ma < 0.1 keV) not accessible to other direct searches. Moreover 
it is a relatively simple technique that can be directly applied to detectors searching for 
WIMPs. 

Primakoff conversion using a crystal lattice has already been employed in two germa- 
nium experiments: SOLAX and COSME-II [^] with the ensuing limits for axion- 
photon coupling gayy < 2.7 X 10-9 GeV-i and gayy < 2.8 x 10"^ GeV"^ respectively. These 
constraints are stronger than that of the Tokyo axion helioscope for nia > 0.26 eV 



and do not rely on astrophysical considerations (i.e. on Red Giants or HB stars dynamics 
||63|| ). The orientation of the crystal was not known so that the data were analyzed taking 
the angle corresponding to the most conservative limit. 

It has been noted that the model that yields the solar axion fiirxes used to calculate 
the expected signals is not compatible with the constraints coming from helioseismology if 
Qa-yy > 10^9 GeV~^ 1^^. This would imply a possible inconsistency for solar axion limits 
above that value, and sets a minimal goal for the sensitivity of future experiments. 

The use of CUORE to search for solar axions via Bragg scattering should have a priori 
some advantages with respect to germanium detectors, because of the larger mass and the 
known orientation of the crystals. On the other hand, as the cross-section for Primakoff 
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conversion depends on the square of the atomic number, Te02 will be a priori a better 
candidate than Germanium. Needless to say that a low energy threshold is mandatory 
because the expected signal lies in the energy region 2 keV < E < 10 keV and is peaked 
at E ~ 4 keV. 

A detailed analysis has been performed ||66| for a Te02 crystal (which has a tetragonal 
structure |^5[) assuming different values for the experimental parameters. As it is shown 
in Ref. the bound on axion-photon coupling which a given experiment can achieve 
can be estimated through the expression: 

9a,, < 9^ - k ( L^^y^) X 10-^ GeV-^ 

^''^^ Vc/keV/kg/dayM T J 

where k depends on the crystal structure and material, as well as on the experimental 
threshold and resolution. For the case of Te02 and a threshold of 5 keV, k has been 
calculated to be = 3.0 assuming an energy resolution of 2 keV. The computation of this 
expression for some assumed values of the experimental parameters is shown in table |^ 
for CUORICINO and CUORE. In all cases flat backgrounds and 2 years of exposure are 
assumed. 

It is worth to notice the faible dependence of the ultimate achievable axion-photon 
coupling bound on the experimental parameters, background and exposure MT: the 1/8 
power dependence of Qa,, on such parameters softens their impact in the final result. The 
best limit shown in table ^ is in fact only one order of magnitude better than the present 
limits of SOLAX and COSME-II. The Qay-f bound that CUORE could provide is depicted 
comparatively to other limits in figure ^. 

The limit which can be expected from the CUORICINO experiment is comparable to 
the helioseismological bound mentioned before. CUORE could go even further (see Figure 
^. It should be stressed that the bounds on ga-y-y obtained with this technique stagnate 
at a few xlO~^° GeV~^, not too far from the goal expected for CUORE, as has been 
demonstrated in [Q. No bounds below these limits can be expected by this technique 
from other crystal detectors like Nal, Ge or Te02, and consequently, there are no realistic 



chances to challenge the limit inferred from HB stars counting in globular clusters ||63 
and a discovery of the axion by CUORE would presumably imply either a systematic 
effect in the stellar-count observations in globular clusters or a substantial change in the 
theoretical models that describe the late-stage evolution of low-metallicity stars. To obtain 
lower values of ga-y-y one should go to the magnet helioscopes like that of Tokio |^ and 



that of CERN (CAST experiment [|6^ currently being mounted). In particular, the best 
current experimental bound of ga,, comes from the Tokyo helioscope: ga-y-y < 6 x 10~^° 
GeV"^ for rua < 0.03 eV and ga-y, < 6.8 - 10.9 x 10"^° GeV"^ for rua ~ 0.05 - 0.27 eV. 
The sensitivity of CAST is supposed to provide a bound ga-^/-y < 5 x 10~^^ GeV~^ or even 
lower. 

Recently the DAMA collaboration has analyzed 53437 kg-day of data of their Nal set 
up |6^, in a search for solar axions, following the techniques developed in ref. ||6^, where 
a calculation of the perspectives of various crystals detectors (including Nal) for solar 
axion searches has been made. The DAMA result gaj-y < 1.7 x 10~^ GeV~^ improves 
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slightly the limits obtained with other crystal detectors ^ and agrees with the result 
predicted in ref. [|6^] . 



6 Conclusions 

We have reported the perspectives of CUORE, a projected massive 760 kg array of 1000 
Te02 bolometers, and of its first stage CUORICINO, with 42 kg of the same crystals 
(currently being mounted), as far as their physics potential to detect various types of rare 
events is concerned. The estimated background and resolution, based on Monte Carlo 
studies, together with the results obtained in recent improvements in the performances 
of the MIBETA experiment and the information obtained from the preliminary tests of 
CUORICINO, have allowed us to assess the potentialities of these experiments for double 
beta decay searches, solar axion detection and WIMP exclusion or identification. In these 
three types of searches, CUORE and to some extent CUORICINO will be powerful tools 
to explore, with higher sensitivity, such rare phenomena. 



7 Acknowledgements 

This work has been partially supported by the Spanish CICYT, contract AEN99-1033 and 
the EU Network Contract ERB-FMRX-CT-98-0167. We are indebted to the CUORE Col- 
laboration for allowing us to use unpublished results of the R&D program of CUORICINO 
and to E. Fiorini for his contribution to the conception and coordination of the CUORE 
project. Useful discussion on the MC simulation of the background with O. Cremonesi, 
M. Pavan and S. Capelli are acknowledged. We also thank G. Luzon for discussions on 
the cosmogenic activation of the crystals. 



References 

[1] N.E. Booth, B. Cabrera and E. Fiorini, Ann. Rev. Nucl. Part. Sci. 46 (1996) 471. K. 
Pretzl, Nucl. Intrum. And Meth. A 454 (2000) 114. E. Fiorini, Nucl. Phys. B Proc. 
Suppl. 91 (2001) 262. L. Mosca, Review Talk given at the IMPF2002, Feb. 2002, Jaca, 
Huesca (Spain), to appear in Nucl. Phys. B Proc. Suppl. 2003. 

[2] E. Fiorini and T. O. Niinikoski, Nucl. Intrum. And Meth. A 224 (1984) 83. 

[3] For a recent survey of WIMP detection, see for instance A. Morales, Nucl. Phys. B 
Proc. Suppl. 87 (2000) 477 ||astro-ph/9912554|| and Review Talk given at TAUP 2001, 



Sept. 2001, LNGS, | |astro-ph/0ll^5^ , Nucl. Phys. B Proc. Suppl. 110 (2002) 39. 



[4] S. Pirro et ai, Nucl. lustrum. Meth. A 444 (2000) 71. 
[5] E. Fiorini, Phys. Rept. 307 (1998) 309. 



19 



[6] A. Giuliani et al, Nucl. Phys. B Proc. Suppl. 110 (2002) 64. 

[7] R. Abusaidi et al. [CDMS Collaboration], Phys. Rev. Lett. 84 (2000) 5699 |astra 



pli/000247li . D. Abrams et al. [CDMS Collaboration], ||astro-pli/0203500i . 



[8] A. Benoit et al. [EDELWEISS Collaboration], Phys. Lett. B513 (2001) 15 j^ito^ 
I ph/010609l . G. Gerbier, talk given at Neutrino 2002, Munich 30 May 2002. 

[9] M. Altmann et al, Proceedings of the X International Symposium on Lepton and 
Photon Interactions at High Energies, 2001, [|astro-ph/010631^ . J. Jochum et al, Nucl. 
Phys. B. (Proc. Suppl.) 87 (2000) 70. 



[10 
[11 

[12 

[13 
[14 



[15 
[16 

[17; 

[18 
[19 

[2o; 



[21 

[22; 



A.Alessandrello et al, Phys. Lett. B486 (2000) 13. 

A.Alessandrello et al., to be published in Nucl. Phys. of Russian Academy of Science, 
||hep-ex/020103g] . 

R. Bernabei et al. [DAMA Collaboration], Phys. Lett. B 480 (2000) 23; R. Bernabei 
et al. [DAMA Collaboration], Phys. Lett. B 450 (1999) 448; R. Bernabei et al, Phys. 
Lett. B 424 (1998) 195. 



S. Cebrian et al, Nucl. Phys. B Proc. Suppl. 110 (2002) 94, [|hep-ex/0111075| . 

A. Morales et al. [IGEX Collaboration], Phys. Lett. B 489 (2000) 268 jhep^ 



ex/0002053|| ; A. Morales et al. [IGEX Collaboration], Phys. Lett. B 532 (2002) 8 



|hep-ex/011006l|] . 
D. Gonzalez et al, Nucl. Phys. Proc. Suppl. 87 (2000) 278. 
P. Belh et al, II Nuovo Cimento A 101 (1989) 959. 
Y. F. Wang et al, Phys. Rev. D64 (2001) 013012. 

M. Aglietta et al, Proc. of 26th Intern. Cosmic Ray Conf., Salt Lake City (USA), 
August 17-25, 1999, HE 3.1.15, [|hep-ex/9905047| . 

C. Arpesella, Nucl. Phys. B (Proc. Suppl.) 28A (1992) 420. 

For an updated review, see for instance A. Morales, Nucl. Phys. Proc. Suppl. 77 
(1999) 335 ||hep-ph/98095i0| ; H. Ejiri, Nucl. Phys. Proc. Suppl. 91 (2001) 255. E. Fior- 
ini, Nucl. Phys. Proc. Suppl. 91 (2001) 262 and O. Cremonesi, Review Talk given at 
Neutrino 2002, Munich, 30 May 2002. 

Q. R. Ahmad et al [SNO Collaboration], ||nucl-ex/0l06m^ . 

S. Pascoli and S. T. Petcov, [|iep-ph / 0205022| . B. Kayser, talk given at Neutrino 
2002, 30 May 2002. 



[23] H. V. Klapdor-Kleingrothaus et al, Eur. Phys. J. A 12 (2001) 147. 



20 



[24 
[25; 

[26; 
[27; 
[28; 

[29 
[30 
[31 

[32; 
[33; 

[34 

[35; 

[36 
[37 

[38; 

[39 
[40 



[41 

[42; 

[43; 
[44; 



X. Sarazin et al, [ |hep-ex/0006031| . 

H. Ejiri et al, Phys. Rev. C 63 (2001) 065501. 

C. E. Aalseth et al, [ ^ep-ex/020102T| . 

H. Ejiri et al, Phys. Rev. Lett. 85 (2000) 2917. 

M. Danilov et al, Phys. Lett. B 480 (2000) 12. 

A. Staudt et a/., Europhys. Lett. 13 (1990) 31. 

T. Tomoda, Rep. Prog. Phys. 54 (1991) 53. 

P. Vogel et al, Phys. Rev. Lett 57 (1986) 3148; Phys. Rev C37 1988 73; M. Moe 
and P. Vogel, Ann. Rev. Nucl. Part. Sci. 44 (1994) 247 

W. C. Haxton et al, Prog. Part. Nucl. Phys. 12 (1984) 409; Nucl. Phys B 31 (Proc. 
Suppl.) (1993) 82; Phys. Rev. D26 (1982) 1085. 

J. Engel and P. Vogel, Phys. Rev. Lett. B225 (1989) 5 

O. Civitarese, A. Faessler, T. Tomoda, Phys. Lett. B194 (1987) 11; T. Tomoda, A. 
Faessler, Phys. Lett. B199 (1987) 473; J. Suhonen and O. Civitarese, Phys. Rev. C49 

(1994) 3055. 

G. Pantis et al, Phys. Rev. C53 (1996) 695 

A. Staudt, T. T. S. Kuo and H. Klapdor, Phys. Rev. C46 (1992) 871 

E. Caurier et al, Phys. Rev. Lett. 77 (1996) 54; Retamosa et al, Phys. Rev. 51 

(1995) 371; A. Poves et al, Phys. Lett. B361 (1995) 1 

J. G. Hirsh et al, Nucl. Phys. A589 (1995) 445; C. R. Ching et al, Phys. Rev. C40 
(1984) 304; X. R. Wu et al, Phys. Lett. B272 (1991) 169; X. R. Wu et al, Phys. Lett 
B276 (1992) 274 

N. A. Bahcall, J. P. Ostriker, S. Perlmutter, P J. Steinhardt, Science 284 (99) 1481. 
G. Jungman, M. Kamionkowski and K. Griest, Phys. Rept. 267 (1996) 195 ||hep- 



ph/9506380 | 



J. Elhs, K. Ferstl and K. A. Olive, Phys. Rev. D 63 (2001) 065016, ||hep-ph/0007TT3 



A. Bottino, F. Donato, N. Fornengo and S. Scopel, Phys. Rev. D 63 (2001) 125003 
||hep-ph/0010203|] . 



L.Bergstrom, Rept. Prog. Phys. 63 (2000) 793, ||hep-ph/0002T26^ . 



A. K. Drukier, K. Freese and D. N. Spergel, Phys. Rev. D 33 (1986) 3495. 



21 



[45] M.L. Sarsa, A. Morales, J. Morales, E. Garcia, A. Ortiz de Solorzano, J. Puimedon, 

C. Saenz, A. Salinas, J. A. Villar, Phys. Lett. B386 (1996) 458. 

[46] D. Abriola, F.T. Avignone, III, R.L. Brodzinski, J.I. Collar, D.E. Di Gregorio, H.A. 
Farach, E. Garcia, A.O. Gattone, C.K. Guerard, F. Hasenbalg, H. Huck, H.S. Miley, 
A. Morales, J. Morales, A. Ortiz de Solorzano, J. Puimedon, J.H. Reeves, A. Salinas, 
M.L. Sarsa, J.A. Villar, Astrop. Phys. 10 (1999) 133. 

[47] P. Belli, R. Bernabei, W. Di Nicolantonio, V. Landoni, F. Montecchia, A. Incicchitti, 

D. Prosperi, C. Bacci, O. Besida, C.J. Dai, Nuovo Czm. C19 (1996) 537. 



[48] F. Donato, N. Fornengo and S. Scopel, Astropart. Phys. 9 (1998) 247 jhep^ 



ph/ 98032951 ; M. Kamionkowski and A. Kinkhabwala, Phys. Rev. D 57 (1998) 3256 
" |hep-ph/9710337|] ; J. D. Vergado s, Phys. Rev. Lett . 83 (1999) 3597; J. D. Vergados, 
Phys. Rev. D62 (2000) 023519 ||astro-ph/000119q] ; P. Ullio and M. Kamionkowski, 
JHEP 0103 (2001) 049 [|hep-ph/0006183|| ; A. M. Green, Phys.Rev. D63 (2001) 043005 
astro-ph/00083Tl] ; C. J. Copi, J. Heo and L. M. Krauss, Phys.Lett. B461 (1999) 43; 



P. Belh et al, [|hep-ph/0203242|| . 

[49] M. Brhhk and L. Roszkowski, Phys. Lett. B464 (1999); P. Belli, R. Bernabei, A. Bot- 
tino, F. Donato, N. Fornengo, D. Prosperi y S. Scopel, Phys. Rev. D61 (2000) 023512 
||hep-ph/990350l|] ; A. Bottino, F. Donato, N. Fornengo and S. Scopel, Phys. Rev. D62 
(2000) 056006 [ |hep-ph/0001309i . 



[5o; 

[51 

[52; 
[53; 

[54 
[55; 

[56; 

[57 

[58; 
[59; 



J. Engel, Phys. Lett. B264 (1991) 114 

Y. Ramachers, M. Hirsch and H. V. Klapdor-Kleingrothaus, Eur. Phys. J. A 3 (1998) 
93. 



F. Hasenbalg, Astropart. Phys. 9 (1998) 339 ||astro-ph/9806T98 



S. Cebrian et al, Astropart. Phys. 14 (2001) 339 ||hep-ph/9912"39l . 

R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38 (1977) 1440. 

K. van Bibber, P. M. Mclntyre, D. E. Morris and G. G. Raffelt, Phys. Rev. D 39 
(1989) 2089. 

R. J. Creswick, F. T. Avignone, H. A. Farach, J. I. Collar, A. O. Gattone, S. Nussinov 
and K. Zioutas, Phys. Lett. B427 (1998) 235 ||hep-ph/9708"2T0l . 



C. Csaki, N. Kaloper and J. Terning, Phys. Rev. Lett. 88 (2002) 161302. 
E. Masso and R. Toldra, Phys. Rev. D 52 (1995) 1755 ||hep-ph/9503"293| . 
E. A. Paschos and K. Zioutas, Phys. Lett. B323 (1994) 367. 



22 



[60] F. T. Avignone, D. Abriola, R. L. Brodzinski, J. I. Collar, R. J. Creswick, D. E. Di- 
Gregorio, H. A. Faracli, A. O. Gattone, C. K. Guerard, F. Hasenbalg, H. Huck, 
H. S. Miley, A. Morales, J. Morales, S. Nussinov, A. Ortiz de Solorzano, J. H. Reeves, 



J. A. Villar, and K. Zioutas, Pliys. Rev. Lett. 81 (1998) 5068 ||astro-pli/ 9708008 



[61] I. G. Irastorza, S. Cebrian, E. Garcia, D. Gonzalez, A. Morales, J. Morales, A. Or- 
tiz de Solorzano, A. Peruzzi, J. Puimedon, M.L. Sarsa, S. Scopel and J. A. Villar, 
Nucl. Phys. Proc. Suppl. 87 (2000) 102 ||astro-ph/991249l|] ; A. Morales, F.T. Avi- 
gnone III, R.L. Brodzinski, S. Cebrian, E. Garcia, D. Gonzalez, I.G. Irastorza, H.S. Mi- 
ley, J. Morales, A. Ortiz de Solorzano, J. Puimedon, J.H. Reeves, M.L. Sarsa, S. Scopel, 
J.A. Villar, Astropart. Phys. 16 (2002) 325 llhep-ex/OlOlU??! 

[62] Y. Inoue, T. Namba, S. Moriyama, M. Minowa, Y. Takasu, T. Horiuchi and A. Ya- 
mamoto, | ]astro-ph/0204338|] . 

[63] For an updated review see, for instance, P. Sikivie, Nucl. Phys. Proc. Suppl. 87 
(2000) 41 and G. Raffelt, Nucl. Phys. Proc. Suppl. 77 (1999) 459. [|hep-ph/9806506 l 

[64] H. Schlattl, A. Weiss and G. Raffelt, Astropart. Phys. 10 (1999) 353|Eepj 
I ph/9807476i . 

[65] M. D. Ewbank, P. R. Newman, E. Ehrenfreund and W. A. Harrison, J. Phys. Chem. 
Solids 44 (1983) 1157. 

[66] S. Cebrian, E. Garcia, D. Gonzalez, I.G. Irastorza, A. Morales, J. Morales, A. Ortiz 
de Solorzano, J. Puimedon, A. Salinas, M.L. Sarsa, S. Scopel, J.A. Villar, Astropart. 
Phys. 10 (1999) 397 ||astro-ph/981 13591] . 



[67] K. Zioutas et al, Nucl. lustrum. Meth. A 425 (1999) 482 |^stro-ph/980TT76 
[68] R. Bernabei et al, Phys. Lett. B 515 (2001) 6-12 



23 



Vorti 




130 rp 

It! 




136Yp 




I.IZ X iU 


-13 


O.oo X iU 


-13 


l.lo X iU 


-13 


l.lz X 10 


-13 


A OA X , in 

4.84 X 10 


-13 


1 o^x, in- 

1.87 X 10 


-13 


1.87 X 10- 


-14 


3.96 X 10 


-13 


7.9 X 10- 


14 


1.54 X 10" 


-13 


1.63 X 10 


-12 






1.13 X 10" 


-13 


1.1 X 10- 


12 






1.21 X 10- 


-13 


5.0 X 10- 


13 


1.73 X 10- 


-13 


7.33 X 10- 


-14 


3.0 X 10- 


13 


1.45 X 10- 


-13 


1.42 X 10- 


-14 


1.24 X 10 


-13 


9.3 X 10- 


14 


5.8 X 10- 


13 


3.18 X 10 


-12 






1.5 X 10- 


14 






2.16 X 10- 


-14 


9.5 X 10- 


14 


3.6 X 10- 


13 


6.06 X 10- 


-14 



nuclear model 



QRPA 
QRPA 
QRPA |1 
Weak Coupling SM |3| 
Generalized Seniority |^ 
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QRPA without pn pairing p5[] 



QRPA with pn pairing [35 
|3§ 



Large basis SM 
Operator expansion method 



Table 1: 2/50i/ nuclear merits (y ^) of emitters used in some source=detector calorime- 
ters, according to various nuclear models. 
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Mass 


Resolution 


Threshold 


Background 


g^^ (2 years) 
(GeV-i) 


(kg) 


(keV) 


(keV) 


(c/kg/keV/day) 


42 


2 


5 


0.1 


1.3x10-^^ 


760 


2 


5 


0.01 


6.7x10-^° 



Table 2: Expected limits on the photon-axion coupling for 2 years of exposure of CUORl- 
CINO and CUORE assuming the quoted values for the experimental parameters 
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Figure 1: Exclusion projected for 2 years of CUORICINO assuming a threshold of 5 
keV, a low energy resolution of 1 keV, and low energy background levels of 1 and 0.1 
c/keV/kg/day respectively. The closed curve represents the DAMA region. The dashed 
line corresponds to the current MIBETA result. 
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Figure 2: Exclusion projected for 1 year of CUORE assuming a threshold of 5 keV, a low 
energy resolution of 1 keV, and low energy background levels of 0.1 and 0.01 c/keV/kg/day 
respectively. The closed curve represents the DAMA region. The dashed line corresponds 
to the current MIBETA result. 
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Figure 3: Sensitivity plot in the (m, a) plane for CUORICINO, assuming a threshold of 5 
keV, flat background 6=1 (solid line) and 0.1 c/keV/kg/day (dashed line) and two years 
of exposure (84 kg year). It has been calculated for (5^) = 5.6 (see the text). The closed 
contour represents the 3cr CL region singled out by the modulation analysis performed by 
the DAMA experiment 



121, and the cross indicates the minimum of the likelihood found 



by the same authors. 
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Figure 4: The solid lines represent the sensitivity plot in the (m, a) plane for CUORE, 
assuming a threshold of 5 keV, two years of exposure (1500 kg year) and flat backgrounds 
of 0.1 and 0.01 c/keV/kg/day. It has been calculated for {S"^) = 5.6 (see the text). 
The sensitivity curve has been also calculated for a possible threshold of 2 keV with a 
background of 0.01 c/keV/kg/day (dashed line). The closed contour represents the 3a 
CL region singled out by the modulation analysis performed by the DAMA experiment 



12 1, and the cross indicates the minimum of the likelihood found by the same authors. 
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Figure 5: Best bound attainable with CUORE (straight hne labelled "CUORE") com- 
pared with others limits. 
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